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ABSTRACT 

It has been shown based on radio variability arguments that some BALQSOs 
(broad absorption line quasars) are viewed along the polar axis (orthogonal to 
accretion disk) in the recent article of Zhou et ah These arguments are based 
on the brightness temperature, exceeding 10 12 0 K which leads to the well- 
known inverse Compton catastrophe unless the radio jet is relativistic and is 
viewed along its axis. In this letter, we expand the Zhou et al sample of polar 
BALQSOs using their techniques applied to SDSS DR5. In the process, we clarify 
a mistake in their calculation of brightness temperature. The expanded sample 
of high Tb BALQSOS, has an inordinately large fraction of LoBALQSOs (low 
ionization BALQSOs). We consider this an important clue to understanding 
the nature of the polar BALQSOs. This is expected in the polar BALQSO 
analytical/numerical models of Punsly that predicted that LoBALQSOs occur 
when the line of sight is very close to the polar axis, where the outflow density 
is the highest. 

Subject headings: (galaxies:) quasars: absorption lines — galaxies: jets — (galax- 
ies:) quasars: general — accretion, accretion disks — black hole physics 


1. Introduction 

About 15% - 20% of quasars show broad UV absorption lines (loosely defined as absorb- 
ing gas that is blue shifted at least 5,000 km/s relative to the QSO rest frame and displaying a 
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spread in velocity of at least 2,000 km/s) outflows (Weymann 1997; Hewett and Foltz 2003; 
Reichard et al. 2003). The implication is that large volumes of gas are expelled as a byprod- 
uct of accretion onto the the central supermassive black hole. Understanding the details of 
this ejection mechanism is a crucial step towards understanding the physics of the central 
engine of quasars. Although evolutionary processes might be related to BAL outflows, it is 
widely believed that all radio quiet quasars have BAL flows, but the designation of a quasar 
as a BALQSO depends on whether the line of sight intersects the solid angle subtended by 
the outflow. The standard model of quasars is one of a hot accretion flow onto a black hole 
and a surrounding torus of molecular gas (Antonucci 1993). The BAL outflow can be an 
equatorial wind driven from the luminous disk that is viewed at low latitudes, just above the 
molecular gas, Murray et al (1995), or a bipolar flow launched from the inner regions of the 
accretion flow (Punsly 1999a,b). The most general feature of 3-D simulations of accretion 
flows onto black holes is huge ejections of gas from the central vortex of the accretion flow 
along the polar axis which rivals the accreted mass flux for rapidly spinning black holes (De 
Villiers et al 2005; Hawley and Krolik 2006). Thus, it is of profound theoretical interest to 
look for evidence of these endemic polar ejecta. 

BALQSOs are so distant that direct imaging of the BAL region is beyond the resolution 
of current optical telescopes. Thus, much of the discussion of BAL geometry is based on 
deductive reasoning. A novel idea for finding the orientation of BALQSOs was developed 
in Zhou et al (2006). They used radio variability information to bound the size of the the 
radio emitting gas then deduced that the radio emission must be viewed close to the polar 
axis and emanate from a relativistic jet, thereby avoiding the well known inverse Compton 
catastrophe. They found BALQSOs that satisfy these conditions in SDSS DR3. We expand 
their methods to SDSS DR5 and apply these results to the physics of the BAL wind launching 
mechanism. 


2. The Brightness Temperature 

Equation (5) of Zhou et al (2006) is not well defined since it is unclear what frame of 
reference is used to evaluate the quantities contained within. Apparently, this gave rise to 
estimates that are (1 + z) 2 3 larger than what we find in equation (2.2), below. In this section, 
we derive a formula for T b . Physically, it is T b evaluated in the cosmological frame of reference 
of the quasar, (T b ) q , that is relevant for assessing the ’’inverse Compton catastrophe.” We 
want to express this in terms of observable quantities at earth designated by the subscript 
”o.” Consider a source in which the monochromatic intensity has increased by an amount 
I{v) q in a time A t q . The brightness temperature associated with the change in luminosity 
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is {T b )g 
(1978) 


(7 (v) q c 2 ) / (2k h v q ) . From the monochromatic version of Liouville’s theorem, Gunn 



(l + z) : 




( 2 - 1 ) 


where {F v ) 0 is the flux density observed at earth. The solid angle subtended by the source, 
Q 0 , is bounded by the causality requirement that the source could not have expanded more 
rapidly than the speed of light during a time, A t q , Q, 0 < (or At q ) 2 / (4(1^), where the angular 
diameter distance is 3a (Gunn 1978). In a a cosmology with i7 0 =70 km/s/Mpc, = 0.7 
and fi m — 0.3, we can use the expression for (1a given by Pen (1999), that is accurate to 
< 1% relative error, along with (2.1) to find: 


(T b ) 


8.0 x 10 12 (1 + z) 


-Z 2 (AiT(mJy)) 0 0 K , 


(2-2a) 


9 (p 0 /lGHz) 2 (A£ 0 /lyr) 2 ' 

Z = 3.31 - (3.65) 

x ([(1 + zf - 0.203(1 + zf + 0.749(1 + zf + 0.444(1 + z) + 0.205] ~°' 125 ) {2- 2b) 


where AFj,(mJy) is the change in flux density in mJy measured at earth at frequency v Q 
during the time interval A t 0 . 


When (T b ) q > 10 12 0 K, the inverse Compton catastrophe occurs (Kellermann & Pauliny- 
Toth 1969). For an unresolved source, (F„) 0 is Doppler enhanced relative to the intrinsic 
flux density, {F u ) 0 = T !+a ( 7b) intrinsic, where the Doppler factor, <5 is given in terms of T, the 
Lorentz factor of the outflow; (3, the three velocity of the outflow; and a, is the spectral 
index, <5 = l/[r(l — (3 cos 6)] (Lind and Blandford 1985). Thus, Doppler beaming can cause 
(Tb) q to be enhanced by a factor of S 3+a in (2.2). Assuming a flat spectral index for the 
unresolved relativistic jet, we choose a — 0. This allows us to find a minimum Doppler factor 
that avoids the inverse Compton catastrophe, S min [(Tb) q ] = [(T&) g / (10 12 0 K)] 0 333 . If the jet 
plasma propagates with S > d m i n [(Tb) g )} then the intrinsic T b in the frame of reference of the 
jet plasma is sufficiently low that there is no inverse Compton catastrophe. 


3. Polar BALQSOs 

Table 1 is a list of BALQSOs in which the radio variability requires that the jet is 
propagating well within 35 0 to the line of sight in order for the jet plasma to satisfy (T b ) q < 
10 12 0 K. Since rather small variations of flux density create these conditions we must first 
prove that the sources are truly variable. We choose the condition for variability to be 

Afp — Smp _ . . 

cw = -=== > +3 , (3-1) 

y'&P'P + <J NP 
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Table 1: Properties of the Radio Variable BAL Quasars 8 


Name 
(SDSS J) 

z 

Sedate 

(mJy) 

Sedate 

(mJy) 

a 

BI CIV 
(km/s) 

bi Aim 

(km/s) 

BI Mgll 
(km/s) 

Type 

081102.91+500724.5 b 

1.84 

23.07: 05/23/97 

18.80: 11/15/93 

5.7 

617 

... 

25 

LoBAL 

081618.99+482328.4 

3.57 

72.29: 05/01/97 

61.10: 11/15/93 

5.2 

2576 

... 

.... 

HiBAL 

081839.00+313100.1 

2.37 

8.82: 10/23/95 

7.2: 12/15/93 

3.1 

? 

... 

.... 

HiBAL 

082817,25+371853. 7 C 

1.35 

21.18: 07/23/94 

14.50: 12/15/93 

10.5 


428 

1564 

LoBAL 

093348.37+313335.2 

2.60 

18.35: 10/23/95 

15.90: 12/15/93 

3.7 

4107 

1405 

.... 

LoBAL 

104106.05+144417.4 

3.01 

27.46: 12/99 

19.00: 12/06/93 

11.4 

2527 

... 


HiBAL 

113445.83+431858.0 

2.18 

27.38: 2/20/97 

24.90: 11/15/93 

3.4 

10643 


1149 

LoBAL 

134652.72+392411.8 

2.47 

3.60: 08/19/94 

2.20: 04/16/95 

3.0 

1505 

557 

.... 

LoBAL 

142610.59+441124.0 

2.68 

6.78: 03/27/97 

5.0: 03/12/95 

3.6 

7 

? 

.... 

LoBAL 

145926.33+493136.8 d 

2.37 

5.22: 04/17/97 

3.60: 03/12/95 

3.4 

9039 

329 

.... 

LoBAL 

155633.77+351757.3 

1.49 

30.92: 07/03/94 

26.90: 04/16/95 

4.3 

e 

... 


FeLoBAL 

165543.24+394519.9 

1.75 

10.15: 08/19/94 

8.50: 04/16/95 

3.0 

? 

... 


HiBAL 


“The first column in table 1 is the source name, followed by the redshift. Columns (3) and (4) are the FIRST 
and NVSS peak flux densities at 1.4 GHz, respectively. Column (5) is the statistical significance of the radio 
variability computed using (3.1). Columns (6) - (8) are the BALnicity indices, Weymann et al. (1991), for 
CIV, A1II and Mgll, respectively. The last column denotes whether the source is a low ionization or a high 
ionization BALQSO. 

6 There is a new spectrum in SDSS with much higher S/N than the spectrum in Zhou et al (2006) that clearly 
shows broad Al III absorption 

c from Zhou et al (2006), however there is now a much higher signal to noise spectrum in SDSS than the one 
used in Zhou et al (2006) and we have recomputed the BALnicity indices 

d from Zhou et al (2006) 

e identified as FeLoBALQSO in Becker et al. (1997) 
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where Spp and S NP denote the peak flux density at 1.4 GHz measured by the FIRST and 
NVSS surveys, respectively, a F p and a F p are the FIRST and NVSS peak flux uncertainties, 
respectively. The beam for the NVSS survey is considerably larger than FIRST ( FWHM ~ 
45" versus FWHM ~ 5”). Thus, a compact source <C 5" (like a compact radio core which is 
the putative site of variable activity) should be detected with equal sensitivity if it exceeds 
the flux limit of the samples. The NVSS beam should pick up more extended flux than 
FIRST. Variable BALQSO fields in NVSS with confusion from nearby sources were omitted 
from table 1. Thus, a FIRST peak flux density larger than an NVSS peak flux density by 3 a 
(hence the plus sign on the RFIS of eqn. (3.1)) is true variability and not an artifact 
of the different beam sizes. One improvement from Zhou et al (2006) is that we used the 
peak NVSS flux density as opposed to the integrated NVSS flux density. This is preferred for 
the comparison described in (3.1), since the variable radio flux is most likely coming from 
sub-arcsecond regions of the jet and this choice is less sensitive to low surface brightness 
noise that can be picked up in the integrated flux (Condon 1997). 

It is crucial to assess the significance of the 3a results in table 1. In order to clarify the 
statistics, we must delineate the super sample from which the sources in table 1 are drawn. 
The potential variable BALQSOs in DR5 that can be detected by equation (3.1) require two 
conditions. First of all, the sources must have a BALnicity index > 0 as defined in Weymann 
et al. (1991) for at least one line either, CIV, A1III, or Mgll. Secondly, S F p > 3.55 mJy. 
The NVSS survey has a peak flux limit of 2.5 mJy in general, but occasionally as low as 2.2 
mJy and a minimum measurement uncertainty of « 0.45 mJy Condon et al (1998). Thus, 
Sp P < 3.55 mJy could never satisfy eqn. (3.1). There are 116 BALQSOs that satisfy these 
criteria and they form the master sample. Eqn. (3.1) is a one sided probability, so if the 
errors are Gaussian distributed as in Condon (1997); Condon et al (1998) then there is less 
than a 20% chance of expecting a variable BALQSO from this sample of 116. However, we 
found 21 variable BALQSOs (12 of which made the table because T B was sufficiently high to 
restrict the line of sight to < 35°). Thus, based on Gaussian statistics our findings are very 
significant. However, if the errors in the tail of the distribution are not Gaussian, but arise 
from unforseen sources of measurement error, then the analysis of Condon (1997); Condon 
et al (1998) does not apply and we have no means to assess the significance of our results. 
Note that there are 4 sources in table 1 with a variability that is significant above the 5a 
level. 

If there is no existing spectrum in the literature, figure 1 shows the SDSS spectra of 
the polar BALQSOs in the regions of the BALs. The SDSS spectra of radio-variable BALs 
were retrieved from the SDSS database and were analyzed using the I RAF 3 software. First, 


3 IRAF is the Image Reduction and Analysis Facility, written and supported by the IRAF programming 
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the spectrum was de-reddened using the Galactic extinction curve, Schlegel et al. (1998), 
then the wavelength scale was transformed from the observed to the source frame. The 
spectra were fitted in XSPEC with a powerlaw plus multiple gaussian model, including 
fits for both emission lines and absorption lines (Arnaud 1996). All the model parameters 
were kept free. The best fit to the SDSS data was determined using x 2 minimization. If 
an emission bump around 2500 A is present then the continuum fit in the CIV (1550 A) 
region was extrapolated to longer wavelength. The BALnicity indices quoted in table 1 
are a consequence of the method of spectral fitting described above and other methods 
might produce different results. However, the exact BALnicity index is not critical to this 
discussion, the sources in figure 1 clearly show BALs and that is the essential point of 
relevance here. 


In table 2, we calculate (' T ^) q , using equation (2.4) and the peak fluxes from table 1. 
From this, we calculate 5 min [(T h )q\ in column (3) as described at the end of section 2. For 
each value of d min , one can vary (3 in the definition of 5 to find the maximum value of 9, 
0max{5min[{T b )q]}, that is compatible with 5 min : 


$max {5min[(T b ) q \} = MaX | 0 


arccos 


( VEIL ] 

(-wKjyj ) 



(3-2) 


The value of 9 max in column (4) is an extreme upper limit because it is the third in a chain 
of bounds. 


• The value of (Tb) 9 in table 1 is a lower bound, since NVSS might pick up some extended 
flux that is missed by the peak FIRST measurement and the radio core was actually 
weaker at the time of the NVSS measurement than indicated by the peak flux density, 
i.e., (AFk(mJy)) 0 is underestimated in (2.4). Also, (2.3) is an inequality. From the 
discussion at the end of section 2, the larger (T&) g , the larger 5 min 

• the larger 5 min , the lower 6 max by (3.2) and the true value of the jet plasma 5 can be 
much larger than 6 m i n [(T b ) g } 

• The actual line of sight to earth satisfies 6 < 9 max by definition. 

Thus the condition for inclusion into table 2, 9 max < 35°, means that the jet is likely to be 
propagating very close to the pole. 


group at the National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
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Wavelength (A) Wavelength (A) 






Wavelength (A) Wavelength (A) 


Fig. 1. — Galactic extinction corrected source- frame SDSS spectra of 8 radio- variable BALs. 
These spectra were fitted in XSPEC with a power-law plus multiple gaussian (both for 
emission lines and absorption lines) model (Arnaud 1996). If an emission bump around 2500 
A is present then the continuum fit in the CIV (1550 A) region was extrapolated to longer 
wavelength. The best fit was determined using % 2 minimization. The model emission spectra 
are plotted in solid lines in magenta color (less the modeled absorption troughs). 








4. Discussion 


Figure 2 is a small modification to Fig. 9 of Punsly (1999b). In section 4.5 of Punsly 
(1999b), it is discussed how a relativistic jet (in red) can coexist nested inside of the bipolar 
BAL wind. This jet can emit radio flux that is beamed to within a small angle of the polar 
axis. Within the axisymmetric version of model, LoBALQSOs exist for lines of sight within 
15° of the polar axis. The density of the wind is highest nearest the polar axis as indicated by 
the gray shading in figure 2 and LoBALQSOs are viewed closer to the polar axis on average 
than HiBALQSOs. The accretion disk coronal X-rays are screened from the BAL wind by 
the dense base of the jet which provides hydrogen column densities of X# ~ 10 25 cm~ 2 . 
Similarly, the different lines of sight near the polar axis, through the inhomogeneous wind, 
naturally provide larger Compton scattering columns and more attenuation of the far UV 
from the inner disk and less attenuation of the near UV and optical from the outer disk, 
making the spectrum appear red. The more extreme and highly polarized LoBALs obviously 
occur in objects without a perfectly axisymmetric distribution of dust as discussed in detail 
in Punsly (1999b). 

Since the LoBALQSOs are viewed closer to the jet axis than the HiBALQSOs in the 
axisymmetric model, geometric arguments imply they should typically have jets with larger 
Doppler factors. Thus, LoBALQSOs should occur in samples of variable BALQSOs at an 
inordinately high rate if the axisymmetric version of the polar model applies to a significant 
subpopulation of the BALQSOs. In our master sample of 116 FIRST, DR5 BALQSOs 
defined in section 3, 69 are HiBALQSOs and 47 are LoBALQSOs. Based on table 1, 4/69 
= 5.80% of HiBALQSOs and 8/47 - 17.02% of the LoBALQSOs have {T b ) q > 5 x 10 12o K 
(this condition is equivalent to 6 < 35°). The sample is small, but the fact that in table 1, 
the LoBALQSO likelihood to have a a large T b is 2.94 times that of the HiBALQSOs tends 
to support the notion the polar BALQSO model represents a significant subpopulation of 
BALQSOs. 

Using radio variability arguments, we expanded on the sample of known polar BALQSOs 
begun by Zhou et al (2006). In the process, we noted that these radio variable BALQSOs 
have an inordinately large LoBALQSOs subpopulation. It is interesting that these properties 
are expected based on an existing detailed theoretical treatment and modeling of bipolar 
BAL winds (Punsly 1999a, b). It would be informative to continue to monitor the FIRST 
BALQSOs at 1.4 GHz, with the matched resolution, in order to find more (T b ) q > 10 12 ° K 
BALQSOs and improve the statistical information in table 2. 
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Table 2: Line of Sight to Radio Variable BAL Quasars 


Name 
(SDSS J) 

T b 

(10 12 ° K) 

^min 

^max 

(degrees) 

Type 

081102.91+500724.5 

5.4 

1.7 

34.9 

LoBAL 

081618.99+482328.4 

39.6 

3.4 

17.1 

HiBAL 

081839.00+313100.1 

11.8 

2.3 

26.1 

HiBAL 

082817.25+371853.7 

125.4 

5.0 

11.5 

LoBAL 

093348.37+313335.2 

20.2 

2.7 

21.6 

LoBAL 

104106.05+144417.4 

8.5 

2.0 

29.4 

HiBAL 

113445.83+431858.0 

5.9 

1.9 

33.7 

LoBAL 

134652.72+392411.8 

73.9 

4.2 

13.8 

LoBAL 

142610.59+441124.0 

11.5 

2.3 

26.4 

LoBAL 

145926.33+493136.8 

8.8 

2.1 

29.0 

LoBAL 

155633.77+351757.3 

66.5 

4.0 

14.3 

FeLoBAL 

165543.24+394519.9 

47.5 

3.6 

16.1 

HiBAL 




BALQSO 


UV and X-rays 
from inner disk 


UV from outer disk 

I H Sr10 23 cm’ 2 


BEL cloud 


/ Dusty Torus 


Dusty Torus 


Fig. 2.— The bipolar wind model of Punsly (1999a, b). The gray clouds of the BAL wind 
has an opening angle of ~ 25° - 30°. Nested inside is a relativistic jet (in red). The highest 
BAL wind densities are near the polar axis. Thus, lines of sight near the polar axis, have 
the maximum attenuation of X-ray and ionizing UV radiation from the accretion flow. This 
implies that lines of sight close to the polar axis are more likely to represent LoBALQSOs. 



